Quantum creep in layered antiferromagnetic 
superconductor 



Tomasz Krzysztori 

Institute of Low Temperature and Structure Research, Polish Academy of Sciences, 
50-950 Wroclaw, P.O.Box 1410, Poland 

E-mail: T . KrzysztonOint .pan. .wroc .pi 

Abstract. In the mixed state of layered superconductor the antiferromagnetic order 
of magnetic ions could possibly create the spin-flop domains along the phase cores of 
the Josephson vortices. We discuss how this feature affects the activation energy in 
the macroscopic quantum tunnelling of the Josephson vortex. We show that the action 
and hence activation energy is rendered temperature dependent so that the quantum 
tunnelling rate becomes temperature dependent below the crossover temperature. We 
show that it in such systems in constant temperature it is possible to occur of two 
regimes of the flux creep, thermal or quantum, dependent on the direction or intensity 
of the applied magnetic field. 

PACS numbers: 74.25.Ha, 74.25Qt, 74.20 We 
1. Introduction 

The discoveries of ternary Rare Earth (RE) Chevrel Phases REMoeSs and RERI14B4 
p] compounds with regular distribution of localized magnetic moments of RE 
atoms have proved conclusively the coexistence of various types of magnetism with 
superconductivity. Intensive experimental and theoretical research has shown that 4f 
electrons of RE atoms responsible for magnetism and 4d electrons of molybdenum 
chalcogenide or rhodium boride clusters responsible for superconductivity are spatially 
separated and therefore their interaction is weak. In many of these systems 
superconductivity coexists rather easily with antiferromagnetic order, where usually 
the Neel temperature is lower than the critical temperature for superconductivity 
T c . Quite often the ferromagnetic order is transformed into a spiral or domain-like 
structure, depending on the type and strength of magnetic anisotropy in the system 
[21 E]. For almost two decades the problem of the interaction between magnetism 
and superconductivity has been overshadowed by high temperature superconductivity 
(HTS) found in copper oxides. However, the discovery of the magnetic order in 
Ru-based superconductors [U El El E] inspired a return to the so-called coexistence 
phenomenon |TJ. The interplay between magnetism and superconductivity was studied 
in d-electron UGe2[8] and ZrZn2[9], where itinerant ferromagnetism may coexist with 
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superconductivity, and in heavy fermion UPd2Al 3 [Tl], where magnetic excitons are 
present in the superconducting phase. The recent discovery of the iron pnictide 
superconductors [15] have trigggered broad interest in the mechanism of the coexistence 
of magnetism and superconductivity in this new class of superconductors [16J. They 
exhibit qualitative similarity to cuprates in that superconductivity occurs upon carrier 
doping (electrons, in this case) of pristine compounds that exhibit magnetism [TT] . 

Among classic magnetic superconductors, the Chevrel phases have been studied 
most intensively. These compounds are mainly polycrystalline materials. However, 
some specific features can be measured only on single crystals. One such effect is a 
two-step flux penetration process, predicted in Ref.[T8l IT9] and later observed solely 
in the antiferromagnetic superconductor (bet) ErRh 4 B 4 [20]. This very interesting 
phenomenon was recently rediscovered in DyMoeSs [23], even though good quality single 
crystals of the classic antiferromagnetic superconductors have been measured for years. 
The specific feature caused by the long-range antiferromagnetic order in the mixed state 
of superconductor might be the creation of spin-flop (or metamagnetic) domain along 
each vortex core. 

Consider an antiferromagnet with two magnetic sublattices as an example. An 
infinitesimal magnetic field applied perpendicular to the easy axis makes the ground 
antiferromagnetic state unstable against the phase transformation to the canted phase 
(spin-flop). On the contrary, if the magnetic field is applied parallel to the easy axis 
the antiferromagnetic configuration is stable up to the thermodynamic critical field Ht- 
When the field is further increased a canted phase develops in the system. Assume 
that in the antiferromagnetic superconductor the lower critical field fulfils the relation 
H c i < \Ht and that the external field, H c % < H < \Ht, is applied parallel to the easy 
axis. Then the superconducting vortices appear in the ground antiferromagnetic state. 
If the field is increased above \Ht the phase transition to the canted phase originates 
in the vortex core because the field intensity in the core doubles the external one. 
The spatial distribution of the field around the vortex is a decreasing function of the 
distance from its center. Hence the magnetic field intensity in the neighborhood of the 
core is less then H?- Therefore, the rest of the vortex remains in the antiferromagnetic 
configuration. The radius of spin flop domain grows as the field is increased. 

Thus, in the considered model there are two distinct types of vortices. Possible 
candidate of such system might be ErBa 2 Cu 3 7 [10l[TT]. This compound has tetragonal 
unit cell with small orthorombic distortion in the ah plane. The Er ions form 
two sublattices antiferromagnetic structure of magnetic moments laying parallel and 
antiparallel to the a direction. The theoretical structure used in the present calculations 
consists of superconducting layers of thickness d s and isolating ones of thickness d i: 
d = di + d s . In the isolating layers, the magnetic moments are running parallel and 
antiparallel to the a direction (easy axis). The magnetic field aligned parallel to the 
conducting planes makes the vortex lattice accommodate itself to the layer structure so 
that the vortex cores lie in between the superconducting sheets. 

The structure of a vortex lying in the ah plane in a layered superconductor with 
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Figure 1. Single Josephson vortex lying in the ab plane along the i-axis (a-axis). The 
SF domain induced along the phase core is shown in the gray area. 



Josephson coupling between adjacent layers resembles the Abrikosov's one except that 
the order parameter does not vanish anywhere [30] . Instead there exists a region, tj 
along the plane and d perpendicular to it, where the Josephson current j z is of the order 
of the critical current. In this region, which plays the role of the vortex core, the London 
model fails. Away from the core the streamlines of the shielding super currents, which 
also represents contours of constant magnetic field, are elliptical except for the zigzags, 
shown in Fig. ([!]), due to the intervening insulating layers. 



2. Macroscopic quantum tunnelling through intrinsic pinning barriers 

A current density j, flowing along the planes perpendicular to the applied magnetic 
field exerts a Lorentz force on the vortices in the c direction so that intrinsic pinning 
barriers are formed on strongly superconducting layers. The experimental evidence of 
quantum tunnelling is based on the fact that the magnetic moment relaxation rate 
exhibits two types of behavior as a function of temperature. Above a characteristic 
temperature T in the thermal activation regime the decay rate is of the Arrhenius type 
T ~ exp (— Uo/kBT). Below T , the decay rate is assumed a priori to be essentially 
independent of temperature F ~ exp (— S/h) and is interpreted as arising from the 
quantum tunneling of vortices through intrinsic pinning potential [25j EHJ E7] . In the 
following we show a considerable change of tunneling rate and crossover temperature 
due to the SF phase formation in the vortex core. 



2.1. General formulae 

We shall assume that the vortex line is a straight string-like object of effective mass m 
per unit length moving in a metastable intrinsic pinning potential V(u) and exposed to 
continuous deformation u(x,t) in the z direction. The magnetic field is applied in x(d- 
axis) direction. Following Caldeira and Leggett [12] the vortex is coupled to a heat-bath 
reservoir assumed to be representable as a set of harmonic oscillators interacting linearly 
with vortex. In the semi-classical approximation the quantum decay rate is calculated 
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as a saddle-point solution (bounce) of the Euclidean action S for the string 
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Here /3 = {ksT) 1 , 77 is the viscosity coefficient and r denotes imaginary time. The 
pinning potential V(u) consists of intrinsic periodic part and the Lorentz potential: 
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where j c denotes critical depinning current. For large current, this potential can be 
expanded around the inflection point to give 
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— w :i and w = ^ \ 2^) 2 ma ^ thought as the width of the barrier 
because V(0) = V(w) = 0, and j c is the critical depinning current. The last term 
in Eq.flTJ) is so called Caldeira-Leggett action [T2], which describes ohmic damping 
produced by the coupling to the heat-bath of harmonic oscillators. The line tension e\ 
is different for vortices in two different orientations in the ab plane. In the semiclassical 
approximation the decay rate is given by the value of the action on a classical trajectory 
obtained from the Euler-Lagrange equations of the motion |13j . 
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The classical trajectory Uo(x), for static solution of Eq.(jlj) independent of r , gives 
the activation energy Uo as the static solution of the following equation in the thermal 
regime T > T cr 

d 2 U r 
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Below the crossover temperature T < T cr a new kind of trajectory, periodic in imaginary 
time, develops. Therefore, u(x,r) can be expanded in the Fourier series with Matsubara 
frequencies 
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Substituting this expansion into Eq.(j3]) and linearizing potential around the static 
solution uo(x) one obtains 
d 2 u r 
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Its solution is easily found to be 
Vq = cosh -2 ( 

Substitution Eq. (6) into Eq. (4) results in the following equation 
ld 2 v n 
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Eqj9] has three discrete eigenvalues: — 1,0, | [32]. The negative one determines the 
crossover temperature 
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When the static solution of Eq.([5]) is used the action So arises only from the elastic and 
pinning terms in Eq.((T]) and the dynamic and inertial terms are cancel out. The action 
So = U hf3 is then given by: 

So = 2.7 7 ( £ld y^ '(l-£) 5/ V (13) 

Below T cr the predominant mechanism by which the metastable state decays is quantum- 
mechanical tunneling. The rate at which this takes place is determined by the action for 
the bounce trajectory. To include contributions of effective mass and viscous damping 
in the bounce action one would have to go beyond the harmonic approximation in the 
pinning potential as employed in Eq.([7]) and substitute for u(x, r) the leading terms in 
the Matsubara frequencies. Calculations could be performed in the vicinity of crossover 
temperature in terms of perturbation parameter (1 — T/T cr ). It can be shown that v n 
are of the order (1 — T/T cr ) 2 and the leading-order contribution in perturbation comes 
from v\ ~ cosh~ 3 (. Substituting this term for u in Eq.((T}) one can obtain the total 
action for the bounce trajectory involving contribution of the inertial mass and the 
viscous damping dynamics. 
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2.2. Evaluation of the parameters 

The above calculations apply to both kind of vortices which posses different effective 
masses, line tensions and viscosity coefficients. We shall calculate these parameters as 
the functions of condensation energy accumulated in the vortex cores. For the stationary 
flux flow the viscous force 77 1| is equal to Lorentz force. The electric field generated by 
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the moving vortex is E = B^, so we get E = ^^-j = pj = pNjj-^j where Pat is the 
normal phase resistivity in the ah plane and H c2 is the upper critical field parallel to the 
layers. Finally, 
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where H r = eiK - 2 /® is calculated from the constitutive relation ei = H c \U)n. 

A moving vortex in the magnetic superconductor can transfer energy to the 
magnetic moments by emitting spin waves. This energy transfer gives rise to magnetic 
contribution to the mass and viscosity of the vortex [28] . However the condition for this 
effect to occur is that the velocity of the vortex should exceed the velocity of the spin 
wave which is not the case during the flux creep. 

The effective mass of the vortex can be deduced from the work of Suhl|31|. He 



derived the core contribution to the inertial mass m = |m P ^ HcMo , where m e denotes 
the mass of the electron and €f the Fermi energy, and the electromagnetic contribution 
coming from the energy of the electric field induced by the moving flux. Simple 
estimation shows that this contribution in layered superconductors is 10 -4 of the core 
contribution. Therefore, 

m = £? 

(Pq-k z e F (In k) 

The vortices in two main orientations in the ah plane have different line tension. These 
ones lying parallel to b direction and those laying in the a direction and created in the 
magnetic field less then \Ht have the line tension Sb and these ones lying in the a 
direction but possessing SF domain have the line tension denoted as e a 

e b = e \n— ; e a = — — + — e ln— — ^\m\2 > ( 17 ) 
a 2 128 n(fi H T + M)d 1 \ A ah 

2 

where en = , — , M is the magnetic moment of the SF domain inside the vortex 
and Ht is the thermodynamic critical field of spin-flop domain formation. In order to 
evaluate the above line tensions it is necessary to calculate first Ht and M. Therefore 
a following argumentation is proposed. At low fields, in the vicinity of the lower critical 
field H c i, the intensity of the field in the vortex core is 2if cl [30]. When the external field 
is increased the field intensity in the vortex core increases because of the superposition 
of the fields of the surrounding vortices. The field intensity in the core must reach Ht 
in order to originate a transition to the SF phase. Thus, taking into account only z 
nearest neighbors we can write for the nonunilateral triangular lattice 

p #T = 2potfci + 4ze b (ln^ K (J- ) + 2K | ^-J^ 

= 2e b + o{e b ), (18) 

where c is the lattice constant. Although there are no precise measurements of spin- flop 
transition in the antiferromagnetic high temperature superconductors, we assume that 
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HoHt 40mT. The typical value of 5.5/zs per Er ion per unit cell in ErBa 2 Cu 3 7 gives 
M « 0, 37 T and fi T « 0, 41 T. Since y? #ci = ^6 and d/£ c « 1 we obtain: 

In f ^ 
gg 1 + _36 Vg^^+MgA^/ 

£f> 128 j n f Anb N 

It is possible now to relate the viscosity coefficient Eq. (fT5]) and the mass of the vortex 
Eq.f JT6|) to its line tension Eq. ffT7|) . 

!h = £ ^ (20) 
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With the use of a simple consideration we can estimate the change of j c due to the 
creation of spin-flop domain along the vortex 

2 
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where the Fourier transform of the compression modulus is given by [33] 
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In the a direction, however, we have an additional contribution from the magnetic 
domain dydz ~ 5ipo/SnBT , so we get 



5dB T 

and finally we get j ca ~ 3j c b. 



Jca - Jc + 128A ^ (23) 



3. Motion of the flux in the quantum regime 

Consider now a hollow cylindrical sample of the radius a and the wall thickness g « a 
(this is so called thin wall approximation). The sample is placed in a magnetic field 
directed parallel to the axis of the cylinder and along the superconducting layers. A 
trapped flux in the system is $ = (B in — B ex ), where B in denotes the field inside the 
hole of the sample and B ex outside the sample respectively. The motion of the flux is 
triggered off by an activation process in which segments of the flux line tunnel through an 
intrinsic pinning potential to the neighboring interlayer spacing. By applying Faraday's 
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law we can easily calculate electric field in the sample due to the change of the trapped 
flux 

1 . djj/jc) 



E 



which is equal to the mean electric field associated with the motion of vortices 
E = (p WLd, where L is the length of the sample in the direction of the applied field, 
d interlayer separation and W is the activation probability per unit volume and unit 
time. The probability for quantum tunneling is given by W ~ exp (— S/K). As was 
shown by Tekiel [31] the activation probability of macroscopic quantum excitations is 
proportional to j 3 . Thus we can assume that W = aoj 3 exp (— S/K). Combining all 
above together we get 

Qt = _ r*(o) exp(-S(x) /h) 

Jx(t) X 3 

where x = j/j c , fl = v 9 o a oJc/(/^o7); 7 = is a geometrical factor determined by 
the shape of a sample. Equation (f21~I) was solved numerically for the input parameters 
shown in Tab.fCQ), and the results are shown in the Fig. (jHJ). As we can see the creep is 



-dx 
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d = 10 9 m 
j cb = 10 12 A/m 2 
mt, = 10~ 22 kg/m 

e F = O.leV 



rjb = 10 kg/(s ■ m) 
\ ab = 25 x 10- 9 m 
A c = 125 x 10~ 9 m 
= 150T 



Table 1. Input parameters 
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Figure 2. Left panel: Normalized creep rate as a function of the temperature. Right 
panel: Quantum creep rate (the rate of the decay of the trapped flux $/$o ~ j/jc) as 
a function of the logarithm of time. Dashed line depicts vortices possessing magnetic 
cores. 



slower when vortices possessing magnetic domain occurred in the system. 



4. Discussion of the results 



The main purpose of the present study is to derive expressions for the bounce 
contribution to the action arising from the inertial mass term, elastic and pinning 
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terms, and the damping term. As was explained previously, in antiferromagnetic 
superconductors it is possible that two distinct types of vortices occur in the sample. 
One that is quite similar to the vortices of nonmagnetic superconductor and the other 
possessing a ferromagnetic-like domain around, and inside, the core. What distinguishes 
both types of vortices in the present calculations is their tension Eq. (ll8j) and critical 
depinning current Eq. (|22p . We consider the limit of a large current which is only slightly 
less then critical depining current. In this approximation we numerically calculated the 
crossover temperature from Eq.( fl2|) and ploted it as a function of (FC), fractional current 
j/ 3c 




0.990 0.992 0.994 0.996 0.998 1.000 



j/jc 

Figure 3. Crossover temperature as a function of fractional current j/j c - Dashed line 
depicts vortices possessing magnetic cores. 



As can be seen the difference of crossover temperatures for both types of vortices 
depends on the critical current, and for the value of the FC equal to 0.99 this difference 
reaches almost 5 K, and then vanishes when FC reaches 1. Fig. (Jl]) shows the activation 
energy U calculated from Eq.( fl3|) as a function of FC. The activation energy differs 
about 1.5 meV for 0.99 of FC and goes to when FC reaches 1. 
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Figure 4. Activation energy Uq (in meV) calculated as a function of fractional current 
j/j c - Dashed line depicts vortices possessing magnetic cores. 



Another important effect is that the action and hence activation energy is rendered 
temperature dependent so that the quantum tunneling rate becomes temperature 
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dependent below the crossover temperature. Moreover, including the damping and 
inertial mass contributions reduces the quantum tunneling rate. From the above 
considerations one can draw yet another conclusion. At constant temperature the flux 
creep regime can be altered. If the direction or intensity of the applied magnetic field 
is changed, it is possible to leap from quantum to thermal creep, and vice versa. Let 




Figure 5. Schematic diagram of the leap from thermal to quantum creep, and vice 
versa. 

us see how it can be performed. Let us increase magnetic field intensity, aligned in the 
a direction, to the point marked A. The system is in the thermal creep regime if the 
temperature of the sample is fixed somewhere in the range T crb > T > T cra , as shown 
in Fig.([S]). Then we change the direction of external field from a to b axis. The system 
now leaps over to the point B and finds itself in the quantum creep regime. Doing the 
same operation in the reverse order one enforces the system to crossover from quantum 
to thermal creep. Yet another scenario is possible. When magnetic field is applied along 
a axis and the temperature is kept constant in the interval T cr & > To > T cra , the increase 
of magnetic field intensity above \H T enforces new vortices to occur, and the leap of the 
system from a quantum to thermal creep regime. This scenario could easily be named a 
"creep valve" . It enables to make the creep rate slower or faster simply by manipulation 
of the magnetic field intensity around its value \Hx- 
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